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The turnover of phylogenetic clades across space is a fundamental biodiversity pattern that may depend on 
long-term evolutionary processes, and that has downstream effects on other aspects of diversity including 
species richness and community structure. Limited niche evolution and limited dispersal are two major 
processes causing spatial restriction, and thus turnover, of clades. We studied the determinants of clade 
turnover within the World's richest floristic kingdom, the Neotropics, using the palm family (Arecaceae) as 
a model. We show that continental-scale clade turnover is driven by a combination of limited niche 
evolution — with respect to temperature and soil tolerances — and limited dispersal. These findings are 
consistent with strong dispersal barriers within the Neotropics, and the observation that some palm lineages 
are most diverse in certain biomes or climates. The importance of such deep-time effects suggest that palms 
might be slow to adapt or disperse in response to anthropogenic climate change. 

Most groups of organisms are geographically restricted in their occurrence, sometimes becoming emblem- 
atic for a region, like the lemurs for Madagascar or the palms for the tropics. This turnover of phylo- 
genetic clades in space is a fundamental biodiversity pattern that has long fascinated biogeographers 1 . 
Two major factors seem to restrict the geographic distribution of clades: limited dispersal and limited P-niche 
evolution (P-niche refers to the niche dimensions that determine a species' distribution in space 2 ). By limited 
dispersal we mean the failure of individuals to reach certain environmentally suitable areas, whereas by limited P- 
niche evolution we mean the failure to evolve adaptations to certain environments. The relative importance of 
these two factors is largely unstudied and likely to vary among taxa, regions and biomes 3 . Of note, geographic 
turnover of clades can lead to other diversity patterns, such as variation in species richness or the functional 
composition and diversity of assemblages, if there is variation in species richness or morphology among clades". 
The factors leading to among-clade differences in species richness have received considerable attention 7 , while the 
factors constraining those clades in space — although a crucial piece of the puzzle — have mainly been discussed 
for a few prominent examples 8,9 . 

Tropical America harbours the most diverse flora of the world 10 , and the origin of this wealth of plant species 
has received intense scientific scrutiny 31113 . Scenarios of Neotropical plant evolution involve effects of both 
dispersal barriers and ecological factors, such as biome conservatism 314 . Influential historical events include 
the long-time separation of South America from North and Central America and the subsequent Great 
American Biotic Interchange, marine incursions into the Amazon basin, and the uplift of the Andes 11,15 . These 
geomorphological reorganizations caused strong changes in dispersal opportunities and habitat distribution, with 
abundant dispersal barriers and steep ecological gradients. Here, we provide new insights into the complex 
history of Neotropical plant diversification by analysing the geographic turnover of clades within a major tropical 
plant family, the palms (Arecaceae). 

Palms are a pantropical group largely restricted to megathermal climates 16 and a textbook example of tropical 
niche conservatism 1,9 . Approaching one third of the family's 2,400 species is found in the New World 17 , with all 
but one of native species being endemic. American palms are distributed between 34°N and 34°S (Fig. la), from 
sea level to 3600 m.a.s.l., and occur under a wide variety of environmental conditions 18 . It is known that the broad- 
scale distributions of palm species depend on both environment — mainly climate — and dispersal 19 . In the 
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Figure 1 | a: Distribution, species richness, and biogeographic regions of palms in the Americas. Green shading indicates the species richness of palms in 1 ° 
X 1° grid cells. Stippled lines indicate the biogeographic "palm regions" suggested by A. Henderson 18 : (1) Mexican region, (2) Caribbean region, 
(3) Central American region, (4) Andean region, (5) Amazon region, (6) Central Brazilian region, (7) Atlantic coastal forest, b-e: Examples of palm 
lineages with a spatially restricted distribution, b: Coccothrinax (14 spp., stem age: 10 mya), c: Syagrus (30 spp., 8.5 mya), d: Ceroxylon (11 spp., 11 mya), 
e: Chamaedorea (77 spp., 19 mya) 18,58 . 



Americas, many palm lineages have a spatially restricted distribution 
(e.g. Fig. 1 b-e). Based on the distribution of palm genera, seven 
biogeographic "palm regions" were proposed as putative centres of 
diversification (Fig. la) 18 . Geographic turnover of palm clades is also 
reflected in strong phylogenetic clustering of palm assemblages in 
"botanical countries" (TDWG level 3 units) 20 . A previous study indi- 
cated that the turnover of larger palm clades (subfamilies) is shaped 
by both P-niche evolution and dispersal limitation 14 . The subfamilies 
Arecoideae and Calamoideae were most species-rich in areas with 
high humidity and fertile soils, indicating limited capacity for 
adaptation to other environments 14 . Species diversity within the 
palm subfamilies peaked in different regions after controlling for 
the current environment, suggesting dispersal-limited diversification 
around different areas of origin 14 . However, the numerically dom- 
inant Arecoideae appeared to be close to equilibrium with the current 
environment, showing little evidence for dispersal limitation at the 
clade level 14 . Since this previous study only analysed a few pre- 
defined clades (the subfamilies), it remains unclear to what degree 
American palm lineages in general are spatially restricted by limited 
P-niche evolution or by dispersal. 

Limited dispersal and limited P-niche evolution should lead to 
predictable relationships between phylogenetic turnover (often 
called "phylogenetic beta diversity" 21 , but see ref. 22) and envir- 
onmental and spatial distances among assemblages (Table 1). 
Analysing these relationships has previously provided important 
insights into community assembly, albeit at much smaller scales than 



in the present study 23 25 . Applied at large spatial scales, this approach 
may complement insights from biogeographic analyses 26 and phylo- 
genetic assemblage structure 20 . At regional scales, phylogenetic turn- 
over is often influenced by speciation (mechanisms M 1 and M 2 in 
Table 1), causing closely related species to have strongly disjunct 
distributions (allopatric speciation) or strongly divergent P-niches 
(adaptive speciation along spatio-environmental gradients) 24 . At lar- 
ger (continental) scales this effect is likely less pronounced, and 
phylogenetic turnover is dominated by limited P-niche evolution 
and/or limited dispersal (mechanisms M 3 and M 4 in Table 1). 

The extent to which limited dispersal leads to predictable relation- 
ships between phylogenetic turnover and spatial distance should 
depend strongly on the environmental heterogeneity of the study 
region. Even spatially close locations may be separated by strong 
dispersal barriers (e.g. mountain ranges, seaways), whereas even dis- 
tant locations may be connected by dispersal corridors (i.e. habitat 
similar to the locations themselves). If a correlation between phylo- 
genetic turnover and spatial distance is indeed due to limited dispersal 
(M 3 in Table 1), it should increase when spatial distance is weighted 
with the occurrence of barriers between locations. Measuring barriers 
is not trivial, since any given habitat (e.g. mountain peaks) may con- 
stitute a barrier for some species (e.g. lowland rainforest species), but a 
dispersal corridor for others (e.g. alpine species). We developed a new 
measure of "barrier distance" that accounts for such differences. For 
each pair of locations, we mapped habitat dissimilarity — relative to 
the habitat at the locations themselves — throughout the study region. 



Table 


1 Determinants of phylogenetic turnover among assemblages (see also 21,47 ) 






Mechanism 


Prediction for environmental dissimilarity (d e ) 


Prediction for spatial distance (d s ) 


Mi: 


Allopatric speciation 


No relationship between phylogenetic 


Phylogenetic turnover decreases with d s 






turnover and d e , if controlled for d s 




M 2 : 


Adaptive speciation along 


Phylogenetic turnover decreases with d e 


No relationship between phylogenetic 




spatio-environmental gradients 




turnover and d s , if controlled for d e 


M 3 : 


Dispersal limitation 


No relationship between phylogenetic 


Phylogenetic turnover increases with d s 






turnover and d e , if controlled for d s 




M 4 : 


Limited P-niche evolution 


Phylogenetic turnover increases with d e 


No relationship between phylogenetic 






turnover and d s , if controlled for d e 
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We then quantified the minimum amount of dissimilar habitat that 
needs to be crossed in order to get from one location to the other (i.e. 
the sum of dissimilarity scores of all grid cells crossed). Deviations 
from a straight line were allowed if this led to a lower barrier distance 
value. This captures the occurrence of dispersal corridors that indir- 
ectly connect locations. Of note, this barrier distance is inherently 
correlated with environmental distance between locations, and it 
can only be compared to pure spatial distance if environmental dis- 
tance is controlled. 

Based on previous knowledge of the geographical ecology of palms, 
we derived the following specific hypotheses for our study: i. Clade 
turnover among palm assemblages in the Americas is determined by 
limited dispersal (M 3 ) and limited (3-niche evolution (M 4 ). Speciation 
patterns (JA l and M 2 ) are, if at all, of minor importance at this scale, 
resulting in an overall positive relationship between both environmental 
dissimilarity and spatial distance and phylogenetic turnover, it. The P~ 
niche evolution effect is strongest (based on the numerical dominance 
of Arecoideae, see above). Hi. Climate, particularly temperature 
extremes and seasonality, is most important for this P-niche evolution 
effect, i v. If an observed relationship between phylogenetic turnover and 
spatial distance is caused by dispersal (M 3 in Table 1), this relationship 
should be even stronger when spatial distance is replaced with a distance 
measure that takes dispersal barriers into account. 

Results 

The turnover of clades among palm assemblages in 1 X 1° grid cells 
(in the following: phylogenetic turnover) was significantly positively 
related to (a) spatial distance, (b) a measure combining spatial dis- 
tance and the presence of barriers in between the cells ("barrier 
distance"), and (c) the overall environmental dissimilarity among 
the cells integrating climate, soils, topography, and remotely sensed 
vegetation structure (Fig. 2a-c). All assessments of significance were 
based on a phylogenetic null model (see Methods). Phylogenetic 
turnover was also significantly positively related to environmental 
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Figure 2 | Relationships between phylogenetic turnover (y-axis, in mya) 
and different distance measures, a: geographic (great circle) distance, 
b: barrier distance, an integrated measure of geographic distance and 
dispersal barriers, c: environmental dissimilarity based on 31 
environmental variables describing climate, soil, topography, and 
vegetation structure, d: temperature dissimilarity based on 1 1 temperature 
variables only. All distances are rescaled to a range of [0, 1]. Dashed lines 
are regression lines. 



dissimilarities that included only one aspect of the environment (i.e. 
temperature dissimilarity [Fig. 2d], precipitation dissimilarity, soil 
dissimilarity, topographic dissimilarity, and vegetation dissim- 
ilarity). When using individual environmental variables (e.g. tem- 
perature seasonality, isothermality, etc. - see Supplementary Table 
SI) to compute dissimilarity, 20 out of the 31 single-variable dissim- 
ilarities significantly predicted phylogenetic turnover. All results of 
the aforementioned models can be found in Supplementary Table S 1 . 

The largest amount of variation in phylogenetic turnover (49%) 
was explained by a model that included barrier distance and dissim- 
ilarity in temperature seasonality (Supplementary Table S2). Models 
that included barrier distance consistently outperformed (i.e. 
explained 1.8-7.4% more variation, median 6.6%) the corresponding 
models with simple geographic distance. The fraction of variation 
independently explained by spatial or barrier distance (R p d) was 
always significantly larger than R pd of the null model (P £ 0.001). 
The fractions of variation that were independently explained by 
environmental dissimilarity (R pc ) were generally smaller when bar- 
rier distance was used than when spatial distance was used, with the 
exception of two environmental variables that anyway showed neg- 
ligible R pe values; however, the median difference was small (0.06%). 
The degree to which the environmental variables co-varied with 
spatial/barrier distance was also highly variable; e.g. the fraction of 
variation explained by both environmental and barrier distance ran- 
ged from 0.2% to 25%, depending on the way environmental dissim- 
ilarity was computed. The different environmental variables also 
differed widely in their independent ability to explain variation in 
phylogenetic turnover (R pe of 0-19% with spatial distance, 0-14% 
with barrier distance). Observed R pe was not always significantly 
larger than R pe of the null model (Fig. 3, Table 2). The single-variable 
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Figure 3 | Variation in phylogenetic turnover among palm assemblages 
explained by environmental dissimilarity and spatial/barrier distance. 

Pairs of bars correspond to different measures of environmental 
dissimilarity (Full = all variables PCA, Soil, Topography, Vegetation 
structure, Temperature, and Precipitation) . The upper bar in each pair is 
based on spatial distance and environmental dissimilarity; the lower bar on 
barrier distance and environmental dissimilarity. R p( j: fraction of variation 
explained by spatial/barrier distance only, R 1T1X : fraction of variation 
explained by either spatial/barrier distance or environmental dissimilarity, 
R pe : fraction of variation explained by environmental dissimilarity only. P 
values indicate the probability that the R pe fraction is not larger than the 
R pe fraction of the null model (Ho). Only P-values < 0.05 are shown, 
x-axis: fraction of variation in phylogenetic turnover. 
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dissimilarities of temperature seasonality, temperature annual range, 
and isothermality consistently had the highest independent ability to 
explain variation in phylogenetic turnover (Table 2). Mean temper- 
ature of the coldest quarter and minimum temperature of the coldest 
month also performed better than any other single-variable dissim- 
ilarity. However, the individual effects of these variables cannot be 
separated, because they are highly intercorrelated (Spearman's | p| s 
0.7; see Supplementary Table S3). Temperature dissimilarity also 
explained more variation independently (R pe ) than any other cat- 
egory of variables (i.e. precipitation, soil, topographic, or vegetation 
dissimilarity), and performed even better than the measure including 
all environmental variables (temperature and non-temperature). 
Some non- temperature single-variable dissimilarities also had sig- 
nificant R pe fractions, albeit with low amounts of explained variation. 
These included soil pH and base saturation, total exchangeable bases 
(only with spatial distance), precipitation of the wettest quarter and 
month, annual precipitation (only with spatial distance), and annual 
mean QSCAT, a microwave backscatter measure of vegetation struc- 
ture. Soil as a group also had a significant effect (P = 0.005) when 
analysed together with spatial distance. 

Discussion 

Our expectation (hypothesis i) that both limited |3-niche evolution 
and dispersal limitation would contribute to the spatial turnover of 
palm clades was strongly supported by the results (Fig. 3). Of note, 
the influence of the two factors cannot be fully separated due to 



spatial autocorrelation in the environment; this is evidenced by the 
fractions of the variation partitioning (Fig. 3). However, the 
independent effects of limited niche evolution (R pe ) and dispersal 
limitation (R pt j) are sufficient to conclude that both factors have a 
significant and quantitatively non-negligible effect on the geographic 
turnover of clades. 

Due to dispersal limitation, lineages often diversify within conti- 
nents or islands, but fail to spread to other environmentally suitable 
landmasses. Within continents, dispersal limitation is more difficult 
to detect, because barriers are not always obvious and differ among 
species depending on their niche 27 . Pure spatial distance can serve as 
a crude proxy for the amount of barriers between two locations 28 , but 
this correlation is imperfect. Of note, pure spatial distance will also 
capture "time-for-dispersal" effects 29 , i.e. a limitation that is not 
imposed by barriers, but by the limited speed at which populations 
can spread across a uniformly suitable landscape. This mechanism 
appears to be relevant sometimes even on evolutionary timescales 30 . 
The significant effect of spatial distance on clade turnover in American 
palms (R p d) can thus be interpreted as evidence for dispersal limitation 
on evolutionary timescales caused either by barriers or by time. We 
were able to gain additional insights by devising a complex measure of 
"barrier distance" that aimed to incorporate the species-specific suit- 
ability of areas as dispersal corridors. The fact that barrier distance 
explains more variation of phylogenetic turnover than spatial distance 
shows that the more complex measure better reflects dispersal barriers, 
and confirms that barriers play an important role for the distribution 



Table 2 Amount of variation in phylogenetic turnover explained by environmental dissimilarity based on separate environmental 
variables. Letters indicate variable category: T = Temperature, P = Precipitation, S = Soil, R = Remote sensing, To = Topography. R pe 
is the amount of variation in phylogenetic turnover explained by environmental dissimilarity independent of spatial/barrier distance. The 
environmental variables were ranked according to the size of R pe ("Rank"). The P-value is calculated as n(R pe ° > R pe )/999, where n(R pe ° > 
R pe ) indicates how many randomized phylogenetic turnover matrices produced higher R pe than the non-randomized matrix in the null 
model. Bold variables have a significant R pe at the a = 0.05 level with both spatial and barrier distance 

Spatial distance Barrier distance 
Variable Category R pe Rank P R pe Rank P 


Temperature Seasonality 


T 


0.198 


1 




0.006 


0.144 




1 


0.012 


Temperature Annual Range 


T 


0.166 


2 




0.002 


0.121 




2 


0.005 


Isothermality 


T 


0.142 


3 




0.002 


0.104 




3 


0.007 


Mean Temperature of Coldest Quarter 


T 


0.080 


4 




0.015 


0.051 




5 


0.034 


Min Temperature of Coldest Month 


T 


0.080 


5 




0.005 


0.052 




4 


0.016 


Soil pH 


S 


0.051 


6 




0.001 


0.036 




6 


0.003 


Soil Base Saturation 


S 


0.040 


7 




0 


0.034 




7 


0 


Precipitation of Wettest Quarter 


P 


0.035 


8 




0 


0.022 




9 


0.004 


Mean Temperature of Driest Quarter 


T 


0.035 


9 




0.059 


0.022 




8 


0.108 


Precipitation of Wettest Month 


P 


0.030 


10 




0.002 


0.018 




10 


0.011 


Annual Mean QSCAT 


V 


0.028 


1 1 




0.001 


0.014 




12 


0.042 


Annual Precipitation 


P 


0.027 


12 




0.009 


0.016 




1 1 


0.052 


Annual Mean Temperature 


T 


0.023 


13 




0.051 


0.010 




14 


0.144 


Mean Diurnal Temperature Range 


T 


0.019 


14 




0.029 


0.01 1 




13 


0.089 


Soil Total Exchangeable Bases 


S 


0.012 


15 




0.024 


0.005 




15 


0.164 


Altitudinal Range 


To 


0.008 


16 




0.145 


0.004 




16 


0.273 


Slope 


To 


0.006 


17 




0.221 


0.003 




17 


0.378 


Soil CaC03 


S 


0.006 


18 




0.229 


0 




24 


0.807 


Mean Temperature of Wettest Quarter 


T 


0.005 


19 




0.317 


0 




25 


0.898 


Max Temperature of Warmest Month 


T 


0.004 


20 




0.277 


0 




21 


0.722 


Precipitation of Coldest Quarter 


P 


0.003 


21 




0.401 


0.003 




18 


0.441 


Annual Mean NDVI 


V 


0.003 


22 




0.486 


0 




31 


0.991 


Mean Temperature of Warmest Quarter 


T 


0.003 


23 




0.394 


0 




26 


0.91 


Soil Cation Exchange Capacity 


S 


0.001 


24 




0.316 


0 




30 


0.929 


Precipitation of Driest Quarter 


P 


0.001 


25 




0.61 


0 




27 


0.913 


Precipitation of Driest Month 


P 


0.001 


26 




0.621 


0 




29 


0.944 


Seasonal Variation of NDVI 


V 


0.001 


27 




0.592 


0 




23 


0.729 


Precipitation of Warmest Quarter 


P 


0.001 


28 




0.553 


0.002 




20 


0.381 


Precipitation Seasonality 


P 


0 


29 




0.704 


0.002 




19 


0.393 


Soil Sand Content 


S 


0 


30 




0.556 


0 




22 


0.442 


Seasonal Variation of QSCAT 


V 


0 


31 




0.878 


0 




28 


0.929 
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of palm clades. This clear result is remarkable since barrier distance (i) 
is based on mean environmental conditions within 1X1° grid cells, 
which cannot capture the occurrence of smaller-scale dispersal bar- 
riers or corridors, (ii) ignores differences in niche breadth among 
species, and (iii) assumes that the environment of a given grid cell is 
optimal for all species occurring there. In spite of these limitations, 
barrier distance as calculated here is clearly better suited for modelling 
dispersal limitation than great circle distance. Our findings likely 
reflect the influence of strong large-scale barriers such as seaways 
(e.g. separation between N and S America, isolated islands of the 
Carribean) and the Andes 3133 , but potentially also of less visible envir- 
onmental barriers (e.g. the savannahs separating the Brazilian Atiantic 
Rainforest from the Amazon). 

The geographic turnover of palm clades in the Americas is clearly 
shaped by P-niche evolution. Our finding that assemblages in differ- 
ent environments tended to be phylogenetically distinct, even when 
spatially close and not separated by measurable barriers, and that 
phylogenetic distinctness increases with environmental dissimilarity 
indicates that niche evolution is overall slow compared to the tempo 
of lineage diversification (Table 1: M 4 ). Effects of limited niche evolu- 
tion (sometimes called "phylogenetic niche conservatism", but the 
definition of this term is controversial 34,35 ) on clade distributions 
have been demonstrated previously 36,37 . There is also evidence show- 
ing that limited niche evolution may, via its effect on clade distribu- 
tions, have shaped patterns of community structure 6 or diversity 4,38 
in terrapins, salamanders, and mammals. Our results supported the 
hypothesis iii that evolution of the climate niche, in particular related 
to temperature, had the strongest effect in American palms (Fig. 3, 
Table 2). It is well known that temperature is important in limiting 
species distributions 39 ; this was previously documented in palms 40 . 
Our results are consistent with the fact that palms are inherently ill- 
suited for tolerating extreme temperatures 9,41,42 . Of note, while water 
availability probably limits the geographic ranges of American 
palms 14,43 , the evolution of those range limits does not seem to be a 
limiting factor for clade distributions (Fig. 3). Although scientific 
evidence has been lacking, the distributions of American palm spe- 
cies are thought to be influenced by edaphic conditions 1419 . This was 
confirmed by our results; furthermore, edaphic range limits appeared 
to evolve slowly compared to the speed of clade diversification, 
although the effect was much weaker than for temperature (Fig. 3, 
Table 2). Slow evolution ("conservatism") of edaphic niches has 
also been shown for other plant groups and regions 44 . In contrast, 
vegetation cover and structure either do not limit American palm 
ranges, or these constraints are phylogenetically labile (Fig. 3). 
Vegetation structure characteristics measured by remote sensing 
have been successfully used to model species distributions and 
diversity 45 , but probably influence palm distributions primarily at 
local scales 19 . 

Our results tentatively suggest that dispersal limitation has an 
equal or higher influence on palm clade turnover than niche evolu- 
tion, rejecting our hypothesis ii. Even if the variation in phylogenetic 
turnover that could be due to either of the factors (R m ) is entirely 
attributed to niche evolution, the influence of both factors would be 
about equal (Fig. 3). The dominance of dispersal limitation is also 
supported by the observation that R pt j > R pe except for the two 
models with highest R pc (Supplementary Table S2). However, com- 
parisons of the magnitude of R pe and R p( j should be interpreted with 
caution since these fractions were affected by different error sources. 
The magnitude of R pe is sensitive to the inclusion of relevant envir- 
onmental predictors and their quality; the magnitude of R pc i depends 
on how well we were able to quantify the presence and strength of 
dispersal barriers between locations. Our results concerning the rela- 
tive importance of dispersal and niche evolution rest on the assump- 
tion that these two error sources were of similar size. Of note, 
dispersal limitation on evolutionary timescales often requires some 
degree of limited niche evolution, since barriers are unsuitable areas 



as defined by the niche 27 . If niches were highly evolvable, barriers 
could easily be crossed by adaptation to the environment on the 
barrier. Any signals of dispersal limitation, such as correlations 
between phylogenetic turnover and spatial/barrier distance, are thus 
also evidence of limited niche evolution. However, this effect — 
mediated by dispersal — is different from the direct niche evolution 
effect (Table 1 M 4 ). 

In general, we consider all our estimates conservative. Uncertainty 
was introduced at various points, including the distribution maps, 
the interpolated environmental predictors, the molecular dating pro- 
cedure, and the artificial species level resolution. Relevant envir- 
onmental predictors could also be missing 46 , although the set of 
included variables was extensive. Neither spatial distance nor our 
measure of "barrier distance" are perfect proxies of dispersal limita- 
tion. Any further improvement in the data sources would likely 
increase the amount of variation explained by dispersal limitation 
or niche evolution. When it becomes possible to model biotic inter- 
actions, the influence of other species on the distribution patterns 
and niche estimates can be added to the equation. 

Besides the effects outlined above, the phylogenetic structure of 
assemblages may be influenced by other factors, notably speciation 
patterns 21,24 (Table 1). At regional scales, repeated allopatric spe- 
ciation across a barrier might lead to a negative correlation 
between phylogenetic turnover and spatial or barrier distance 47 . 
This pattern was found in a study on phylogenetic turnover among 
hummingbird communities in Ecuador (ca. 700 km in extent) 24 . 
Since allopatric speciation is the outcome of dispersal limitation, 
the effect of dispersal on the relationship between phylogenetic 
turnover and spatial/barrier distance appears to reverse in dir- 
ection between regional and continental scales. However, the effect 
of allopatric speciation on phylogenetic turnover should also be 
relatively ephemeral 48 ; as soon as sister lineages on each side of the 
barrier start to diversify into clades, phylogenetic turnover 
becomes positively related to barrier distance. Cases of relatively 
recent allopatric speciation are documented in American palms 31 
but clearly do not dominate our results. This might partly be due to 
the lack of terminal (intrageneric) phylogenetic resolution and/or 
our choice of phylogenetic turnover metric, which is known to 
emphasize 'basal' phylogenetic relationships 49 . However, given 
the large extent of the study and the long presence of palms in 
the Americas, it is likely that any signal of allopatric speciation 
would be weak compared to the diversification of clades within 
regions that are delineated by barriers 18 . Still, there might be some 
effect of allopatric speciation that decreases the correlation 
between phylogenetic turnover and spatial/barrier distance. A sim- 
ilar logic applies to adaptive speciation: if sister species segregate 
along spatially structured environmental gradients 50 phylogenetic 
turnover might actually decrease with environmental dissimilarity 
(M 2 in Table 1). Again, this is most likely to occur at spatial scales 
much smaller than the scale of our study. In any case, our inter- 
pretations are conservative: if allopatric or adaptive speciation 
influenced phylogenetic turnover in our study system, it would 
bias the slopes and fractions of variation presented here towards 
zero and insignificance. 

Our findings allow some thoughts on the future of American palm 
diversity under climate change. If niches evolve slowly even on evolu- 
tionary timescales species will not readily adapt to rapidly changing 
environments. Climate-related genetic variation within species 
appears to be widespread 51 , and rapid adaptation to climate change 
occurs 52 . Adaptation at the population level through selection acting 
on existing genotypes may be rapid, but adaptation to truly new 
environments is limited by the supply of new genotypes through 
mutation 51 . The former is conceivable even for groups showing strong 
phylogenetic niche conservatism, but only until the climatic changes 
extend beyond the range of conditions that is covered by existing 
intra-specific niche variation. At the end of the day, adaptation to 
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environments that are not already part of the species' niche may 
be the key to surviving dramatic, rapid climate change. Our data 
suggested that this type of adaptation might be limited in 
American palms; however, more studies are required to improve 
our understanding of the adaptive capacities within this charismatic, 
ecologically and economically important plant group towards a chan- 
ging world. 

In conclusion, we found that the turnover of clades within 
American palms is controlled by a combination of two factors: 1) 
tolerances to temperature extremes and seasonality, and soil prefer- 
ences, evolve slowly compared to the tempo of lineage diversification 
and 2) dispersal is limited on evolutionary timescales, definitely by 
barriers and possibly by time. The results corroborate that phylogen- 
etic history is pivotal for understanding large-scale diversity patterns, 
especially since the turnover of clades can have important down- 
stream effects on patterns of species richness and functional divers- 
ity. Future studies on other groups of organisms should also take into 
account that diversity patterns may be strongly shaped by restricted 
lineage mobility during diversification. Of note, the balance of dis- 
persal limitation and phylogenetic niche conservatism, as well as the 
relevant niche factors, might be very different in other taxa due to 
ecological and life history differences, and in other regions due to 
differences in geomorphology, environment, and history. We found 
that 'barrier distance', our measure of dispersal barriers between 
locations, is a more relevant predictor of phylogenetic turnover than 
spatial distance. This measure should be similarly useful for studies 
on species turnover. Finally, because their dispersal and niche evolu- 
tion are limited on evolutionary timescales, palms might fail to 
migrate or adapt in response to rapidly changing climates. This 



would have serious consequences for tropical biodiversity and for 
cultures that depend on palms as a resource 53,54 . 

Methods 

An overview of the analyses performed for this article is shown in Fig. 4. 

Distribution data. We used distribution maps of all American palms 18 , excluding 
three that are endemic to remote islands with no other palms (Brahea edulis, Juania 
australis, Sabal bermudana). This resulted in a total of 547 species. Due to taxonomic 
developments, close to 800 species of palms are currently recognized for the 
Americas 17 . However, there are no comparable distribution data that match the more 
recent taxonomy, and we are not aware of any bias in taxonomic changes since 1995 
that might influence our results. We registered presence or absence of the included 
547 species within 1° X 1° grid cells {n = 1701 cells with palms, Fig. la). We preferred 
this grid to an equal area grid because it has previously been used in analyses of palm 
species richness 14,55 , turnover 56 , and net diversification 57 . Thus, we accepted small area 
differences among grid cells for the sake of comparability. The richness per grid cell 
ranged from 1 to 83 species (Fig. la) 14 . These cells were used as assemblages in the 
analysis of phylogenetic turnover. 

Phylogenetic tree. We modified a recently published dated version 58 of a genus level 
phylogenetic tree of all palms 59 to fit our assemblage dataset. We deleted genera that 
did not occur in the Americas and replaced terminal branches that represented genera 
with > 1 species with polytomies that included all species in our dataset. The crown 
age of those polytomies had to be chosen arbitrarily; to test the influence of this 
parameter on our results we produced four versions: crown ages equalling 10%, 50%, 
and 90% of the genus stem age, and crown ages computed with a simulation 
approach 60 . The variation caused by to those different assumptions was orders of 
magnitude smaller than the values of interest, and did not influence the qualitative 
outcome of our study. Thus, we present only the results for genus crown ages set to 
50% of the stem age. Other related measures of phylogenetic community structure 
have been found equally insensitive to this parameter 20 . 

Phylogenetic turnover. For all pairwise combinations of grid cells, the comdist 
metric 48 was calculated in R (http://www.r-project.org/). We modified the comdist 
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function of the package picante 61 to improve computational performance {see 
Supplementary Methods online). Comdist is defined as: 

j S s 

comdist ab = — d vPiApjB, ( 1 ) 

i=l ;=1 

where p (j4 is the presence { — 1 ) or absence ( — 0) of species i in assemblage A; d^ is the 
phylogenetic distance between species i and;'; A is the species richness of assemblage 
A. In our case, da was proportional to pairwise divergence time; therefore, comdist 
represented the time in the past when the members of assemblage A diverged, on 
average, from the members of assemblage B. 

Environmental variables. We extracted environmental information for each grid cell 
from Geographic Information System (GIS) layers with Arclnfo 9.3 (ESRI Inc., 
Redlands, California, USA). We obtained a total of 19 climatic variables, 6 soil 
variables, 2 topographic variables, and 4 measures of vegetation structure (see 
Supplementary Table S4 for information on the sources). All variables were averaged 
for each grid cell, except elevation, where a range was calculated to reflect topographic 
heterogeneity. All variables were standardised and centred. 

Environmental dissimilarity. First, we computed a measure of environmental 
dissimilarity where all variables were included on an equal basis. We subjected the 
dataset to a principal components analysis (PCA) to remove multicollinearity 
between variables; all components with standard deviations < 10% of the standard 
deviation of the first component were omitted. Environmental dissimilarity was then 
calculated as the Euclidean distance in the space spanned by the PCA axes. Second, we 
repeated this procedure for five separate sets of variables, including temperature 
(bioclimatic variables 1-11), precipitation {bioclimatic variables 12-19), soil, 
topography, and vegetation structure {QSCAT and NDVI). Third, we computed a 
separate dissimilarity matrix for each environmental factor, with Euclidean distances 
taken directly from the standardised variables. This resulted in a total of 37 matrices of 
environmental dissimilarity for all pairwise combinations of grid cells. 

Dispersal limitation. We computed two different measures to represent the degree to 
which dispersal among grid cells was limited. First, we calculated geodesic distance 
between grid cell midpoints of all pair- wise combinations {'spatial distance'). Second, 
we calculated a 'barrier distance' that takes the environmental insuitability of 
interjacent areas {i.e. dispersal barriers sensu ref. 27) into account. For each grid cell, a 
GIS layer was calculated for the entire study area representing how different the 
environment is from the environment of the focal grid cell ('cost surface'). For each 
pair of grid cells, environment -based dispersal distance was calculated as the average 
length of the least-cost path between them, based on their respective cost surfaces. 
Details on the calculation can be found in the Supplementary Methods online. 

Analyses. We used ordinary least squares linear regression models to determine 
whether environmental dissimilarity (d e ) and spatial/barrier distance (d d ) could 
explain the variation in phylogenetic turnover (d^). As dj we used either spatial 
distance or barrier distance. Three models were built, with d^ ~ d e> d^ — d d , and d^ ~ 
d e + d d . All variables were centred and standardised to obtain comparable effect sizes. 
We used variation partitioning 52 to divide the variance of phylogenetic turnover into 
four fractions, according to whether they could be explained by environmental 
dissimilarity alone (R pe ), spatial/barrier distance alone (R pd )> environmental 
dissimilarity and dispersal distance (R inx ), or unexplained (R un ). These fractions were 
computed from the R 2 values of the d^ -— d e model (R e ), the d^ ~- d d model {R d ), and 
the multiple predictor model (R t ) as follows: R pe = R t — Rj; R pd — R t — R^; R^ = Rg 
+ Rd - R t ; and = 1 - R t . 

Significance testing. Since pairwise measures of dissimilarity are not statistically 
independent, we used a permutation procedure to assess significances. This allowed 
us to account simultaneously for a second problem, namely that phylogenetic 
turnover between assemblages is confounded by the number of shared species (i.e. 
species turnover) as well as differences in species richness. For the purpose of our 
study, we were only interested in the turnover of higher- level lineages, not species. We 
therefore employed a null model where we re-shuffled the species labels on the 
phylogeny 999 times and re-calculated phylogenetic turnover from these randomised 
phylogenies. For each of these null matrices, we performed the regression analysis and 
variation partitioning described above. In these null analyses, comdist still correlated 
with environmental dissimilarity and dispersal distance due to shared species and 
richness differences. The proportion of times that the R pe and R p d fractions, 
respectively, of the null analyses were higher than the R pe /R p a fractions of the real data 
was used as the probability of the hypothesis H 0 that niche evolution/dispersal 
limitation did not affect phylogenetic turnover. This procedure has been previously 
used to test for relationships between phylogenetic turnover and altitudinal difference 
independent of species turnover 23 . 
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